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Abstract 
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1 Introduction 



Even though QCD yields a remarkably good description of the strong interaction, the low 
energy hadron physics has to be modelled phenomenologically. This is due to the fact that 
the usual perturbative approach in the coupling constant cannot be applied to QCD below 
energies of the order of 1 GeV. Most of the phenomenological results were based on PCAC 
and Current Algebra. However, in 1979, Weinberg showed how to reobtain many of these 
predictions by means of an effective lagrangian. 

The fields in that lagrangian are the light mesons, (pions, kaons and etas) which are 
understood as the Goldstone Bosons (GB) arising from the spontaneous breaking of chiral 
symmetry. The lagrangian is built as an expansion in derivatives, that respects the symmetry 
breaking pattern of QCD. Indeed, the first term in the expansion is fixed by the symmetry 
requirements and accounts for the Current Algebra results. The next terms in the expansion 
produce further corrections, which depend on several phenomenological parameters but are 
always consistent with the QCD symmetry constraints. These techniques were later devel- 
oped to one loop in a set of papers by Gasser and Leutwyler 0, ^. They showed how to 
obtain amplitudes involving light mesons, as functions of their momenta, their masses and 
those few phenomenological parameters. 

By fitting these parameters from a few low energy experiments it is then possible to 
obtain successful predictions for other processes. The whole approach is known as Chiral 
Perturbation Theory (ChPT). 

Very recently some partial higher order calculations [Q have appeared in the literature 
as well as a complete two loop calculation of vrvr scattering ^ , which will be needed in order 
to analyze new data to come from DA$NE and Brookhaven. For a general review of the 
available experimental data on pion physics and future prospects, we refer the reader to 0. 

Nevertheless, there are some intrinsic limitations when applying ChPT, namely, the fact 
that the amplitudes calculated within the chiral approach are only unitary in the perturbative 
sense, that is, up to the next order in the external momenta. Such a breakdown of unitarity 
is most severe at high energies, where the external momenta is no longer a good expansion 
parameter, although it can also occur at moderate energies 0. As a result, it is not possible 
to reproduce resonant states, which are one of the most characteristic features of the strongly 
interacting regime. Many different methods have been proposed in order to improve this 
behavior and thus to extend the applicability of ChPT to higher energies; among them: The 
use of Pade approximants 0, the explicit introduction of resonances |T^, the K- matrix 
||1 1|| , the large limit (N being the number of GB) or the inverse amplitude method 

(lAM) 



This work is devoted precisely to the last method, which can be justified within a disper- 
sive approach and can easily reproduce the two lightest resonances: the p(770) in tttt scat- 
tering PI and the K*{892) in irK scattering [|13|. But not only that, the lAM also improves 



considerably the fit to data even in non-resonant channels, almost up to the first two particle 
inelastic threshold (The many particle inelastic thresholds can be neglected since they are 
suppressed by phase space factors). This fit provides a remarkably good parametrization 
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that can be used for other processes. Indeed, in a previous work []T5[, the authors showed 
how it can be used together with a simple unitarization prescription to obtain successful 
results on 77 —>■ n^n^ up to 700 MeV. 

Of course, it is also possible to obtain very good parametrizations 0, |10| of tttt or ttK 
elastic scattering by including all resonant states explicitly. However, our aim choosing 
the lAM is to reproduce these phenomena just with the few phenomenological parameters 
present in the ChPT lagrangian. In this way, even though their masses and widths will not 
be obtained with great accuracy, resonances can be regarded as real predictions. That is one 
of the relevant features of the lAM since other very popular unitarization methods are not 
able to reproduce resonances unless they are explicitly introduced in the calculation. That 
is, for instance, the case with the K-matrix. 

The purpose of this work is, first, to study how high in energies the lAM yields good 
results and what are its limitations. We would also like to know whether it is possible 
to reproduce further resonance states. It is clear that the best candidates are the lightest 
resonances whose dominant decay modes are tttt or ttK. We have listed them in Table 1. In 
case these resonances were not accommodated after our unitarization, it would be interesting 
to understand why. Second, once we have a good fit to these resonances, we want to make 
a complete numerical analysis of several low-energy quantities of interest, like the chiral 
parameters or the scattering lengths, including estimations for their errors. As we will see 
below, we expect that the lAM somehow will include effects that cannot be obtained from 
the pure p"^ expansion. 



Name 


I,J 


Mass 


Width 


Dominant decays 


p(770) 


1,1 


768.8 ± 1.0 


150.3 ± 1.0 


TTTT, 100% 


/o(980) 


0,0 


980 ± 10 


40 to 400 


TTTT, (78.1 ±2.4)% 

KK, (21.9 ± 2.4)% 


/2(1270) 


0,2 


1275 ± 5 


185 ± 20 


TTTT, (84.7 ±2.6)% 


X*(892)± 
K*{892f 


1/2,1 
1/2,1 


891.59 ± 0.24 
896.10 ± 0.28 


49.8 ± 0.8 
50.5 ± 0.6 


ttK, ^ 100% 



Table 1: Lightest resonances with tttt or ttK dominant decay modes. Data taken from |If 



Finally, we would like to comment on another motivation of the present work, which at 
first may not seem very related to the main topic. The philosophy of the chiral approach has 
also reached the description of the strongly interacting symmetry breaking sector (SISBS) 
of the Standard Model . The scalar sector of such a model displays the same symmetry 
breaking pattern as two flavor massless QCD. Hence it is possible to build an effective 
lagrangian, much as it is done for ChPT |[1^ . Although the electroweak GB are not physical, 
using this lagrangian it is possible to obtain predictions for the scattering of longitudinal 
gauge bosons at future colliders, like the LHC. Indeed, there are already experimental 
proposals to measure the electroweak chiral parameters at CMS |^0|. Most of the works on 
the SISBS make use of the Equivalence Theorem |]T^, which allows us to read the observable 
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amplitudes, in terms of longitudinal gauge bosons, directly from those with GB. This theorem 
has been recently proved in the chiral lagrangian formalism and seems to be severely 
constrained by the lack of unitarity. At this point is when the unitarization procedures come 
into play and it is crucial to know whether they are reliable, since what we are now looking 
for are real predictions and not elaborated fits to still unavailable data. 

In Section 2 we review some basic aspects of exact and perturbative unitarity and we 
define the partial waves in elastic scattering. Section 3 introduces the lAM, first with a 
derivation from Dispersion Theory and then by studying the constraints to its applicability. 
Section 4 and Section 5 are organized in the same way, although they refer to SU{2) and 
SU{3) ChPT, respectively: First we apply the JAM to ChPT with the chiral parameters 
obtained from low energy experiments in order to study the I AM predictive power. Next, 
we present an lAM fit to the data. For the best SU{3) fit we present the unitarized results 
for the scattering lengths and some other phenomenological parameters. Then, in Section 6, 
we study the analytic structure on the complex plane of the lAM amplitudes. In Section 7 
we present the conclusions. There is also an Appendix where we give the elastic scattering 
formulae used in this work, as well as a discussion on perturbative unitarity. 



2 Partial waves, phase shifts and unitarity. 

When dealing with strong interactions, it is usual to project the amplitudes in partial waves 
with definite angular momentum J and isospin /, as follows 

tij{s) = f^d{cose)Pj{cose)Tj{s,t) (1) 

where = 2 or 1 depending on whether the particles in the process are identical or not. 
The acceptable isospin values also depend on the process, namely I = 0,1,2 for vrvr elastic 
scattering and I = 1/2,3/2 for ttK. For both reactions the definite isospin amplitudes Tj 
are obtained from a single function. In the first case 

To(s,t,-u) = 3A{s,t,u) + A{t, s,u) + A{u,t, s) 
Ti{s,t,u) = A(t, s,u) — A{u,t, s) 

T2{s,t,u) = A{t,s,u) + A{u,t,s) (2) 

whereas for ttK scattering we can write 

3 1 

Ti/2(S, t, U) = -T3/2(m, t, S) - -T3/2(S, t, u) (3) 

In order to deal with both processes on the same footing, we will label the particles in 
the reaction as a and (3. Thus the Mandelstam variables will satisfy: s + t + u = 2(M^ + M|) 
and the threshold will be at Sth = {M^ + Mp)"^. As it is well known, whenever s > Sth, and 
below inelastic thresholds, the unitarity of the S-matrix imphes 

Imtij = aap \tij\^ (4) 
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where aafs is the two particle phase space. Exphcitly: 



\ 



1 



(A4 - Mfs) 



(5) 



As a consequence of Eq.^ the partial wave can be parametrized as follows: 

1 



tij{s) 



z'^'^^'^ sin57j(s) 



(6) 



and 5ij{s) is called the IJ phase shift. 

We have already mentioned that the ChPT amplitudes are obtained as an expansion in 
external momenta and masses. That is 



(1) , +(2) 



where, for the cases we are interested in, tf] is C(p^), t/j is C(p^), etc... In practice, we can 
only obtain the few first terms of the series above and therefore the amplitude only satisfies 
the unitarity condition perturbatively 



(7) 



Imt 
Imt 



Im(tS + 1 



(0) 

IJ 

(1) 
IJ 







cr. 



(4f + 2t 



I J ^^'^IJ 



CI. 



a(3 



AO) 
^IJ 



.(1) 12 
^IJ I 



The 0{p'^) terms were given by Weinberg and they are called the low energy theorems. 
The next order contributions to tttt scattering were given in 0, H]. The calculation for ttK 



scattering can be found in |2^, |23|, although we have found that the formulae in the literature 
do not satisfy Eq.^. We will comment on that later. Very recently it has appeared the 
complete calculation of the OijP) contribution to elastic vrvr scattering [§. Although we will 
not use it, we will compare some of its results with those of our method. 



3 The Inverse Amplitude Method 
3.1 Derivation from Dispersion Theory 

Let us briefly review the standard derivation [|^, |13[ of the inverse amplitude method, since 
we will use it later in order to understand the applicability of the method. 

Any partial wave obtained from a relativistic Quantum Field Theory should present 
a characteristic analytic structure in the complex s plane. Indeed, the reaction threshold 
becomes a cut in the from sth to +cxo. Due to crossing symmetry, there should be another left 
cut along the negative axis. If we now apply Cauchy's Theorem to our complex amplitudes 
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we obtain integral equations known as dispersion relations. For instance, a three times 
subtracted dispersion relation is 

tiAs) = Co + C,s + C,s' + - r ^ ^^(^^^^ 

vr J{Mc,+M0)-2 s"(s' - s - le) 

Where we have not written explicitly the left cut (LC) contribution. The number of sub- 
tractions needed depends on how the amplitude behaves at infinity in order to ensure the 
vanishing of the contributions coming from closing the integral contour. We have chosen 
three subtractions since we are going to use 0{p'^) ChPT amplitudes which at high s be- 
have as s^. But our arguments remain valid for 0{p^) amplitudes when using four times 
subtracted dispersion relations, etc... 

The ChPT partial waves present both cuts and we can calculate both the subtraction 
constants Co, Ci, C2 and the integrand inside Eq.^ perturbatively 

tfj = ao + ais 

4'j = (,„+M+M^+^ r (10) 

7C J{M^+Mi3)2 s'^{s' - s - te) 

Where we have expanded the subtraction constants in terms of M^/F^. 

The lAM is based on the fact that the function displays the very same analytic 

structure of tjj, apart from some possible pole contributions. For later convenience, we 
will make use of G{s) = tfP /tij. Notice that we have multiplied 1/t by a real function 
without singularities; thus we keep the same analytic structure and we can write a very 
similar dispersion relation: 

2 , r ImG(s')rfs' 



G{s) = Go + Gis + G2s' + - , . ^ + LC{G) + PG (11) 

where PG stands for possible pole contributions. The advantage of using G{s) is that, using 
Eqs.H and ||, we can calculate exactly the integral over the right cut (but not on the left, 
since those equations only hold on the elastic cut), as follows: 

toG = -4Tl^^ = -4°iV = -ImiS (12) 

I ^IJ I 

Note that we are denote by tjj the exact amplitude, which is unknown, although we know its 
analytic properties. In contrast, the expressions for tf] and t^j, etc... have been calculated 
explicitly. 

As we did before, we can also expand the Gi subtraction coefficients in powers of M^/F^, 
and then rewrite the dispersion relation for G{s), which now reads 

^(0)2 

~ Oo + ais — bo — bis — b2S^ 

£!r ■■;f'<^')^^' + (13) 

IT J{Mc+Mfi)2 S'^{s' - S - le) ''-^^ JJ IJ 
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where we have approximated ImG ^ — Imtj^j on the left cut and we have neglected the pole 
contribution. In other words, 

^(0)2 
f-ZJ ''IJ 

This is the JAM result that we are going to use in the present work. Incidentally, Eq.|l4| can 
be understood as the formal [1, 1] Fade approximant of the ChPT amplitude. 



It is important to remark that if we expand again Eq.|14| at low energies, we find 



,(0)2 

tiJ ^ ;(5r4T) ^ + tf] + 0{p') (15) 
tij tjj 

That is, we recover the ChPT expansion. Hence, up to 0{p^) our method and ChPT yield 
the same low energy results if the same chiral lagrangian parameters are used. 

3.2 The applicability of the Inverse Amplitude Method 

Let us review all the approximations made in the previous section, in order to comment how 
they will constraint the lAM applicability: 

3.2.1 The left cut 

In Eq.|13| we have replaced the G{s) left cut integral by that of — t^j (s). As we have remarked 
in the preceeding discussion, Eqs.§ and | are only exact on the right cut. On the left cut we 
cannot write the chain of equalities that lead to Eq.|12|. Nevertheless, if we use the ChPT 
result as an approximation: 

ImG = -tfj'^^ ^ -Imt^^) + Oip') (16) 
I tiJ \ 

we get 

Notice that, in order to obtain the IJ phase shifts, we are going to calculate t/j(s) for real 
s > 4M^. That means that the denominator {s' — s — ie) inside the integrals is never going 
to be very small, which somehow will wash out the error on the left cut. But note also that 
treating differently the right and left cuts violates crossing symmetry. 

Indeed, in it has already been pointed out that the Pade approximants do not 



reproduce correctly the subleading logarithms that would appear at next order in the chiral 
expansion (C(p^) in this case). Of course they would be obtained if we applied the lAM to 
the chiral amplitudes at 0{p^), but still the method would not yield the correct logarithms 
at 0{p^) and so on. At high energies chiral logarithms are not so relevant, but at low 
energies they are a very important feature of ChPT and indeed they can give the dominant 
contribution in some channels. 
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Nevertheless, from Eq.|T5| we see that at low energies the lAM yields the very same 
ChPT expansion, including the dominant chiral logarithms. The contribution from the left 
cut and subleading logarithms is 0{p^). As a consequence, if we try to make a low energy 
fit to the data, the parameters that we would obtain with the lAM would not lie very far 
from those of ChPT, but they will not be the same. That is the reason why, in the following 
sections, we will denote with a hat the parameters obtained from any lAM fit.0 



3.2.2 Resonances and the pole contribution 

In Eq.0, we have neglected the contributions coming from zeros in the amplitude, that will 
appear as poles of the inverse function. There is no way to know a priori whether or not a 
partial wave will vanish for a given value of s, although it is known that chiral amplitudes 
have zeros below threshold, which are known as Adler zeros. Their position is not known 
except for the J = 1, J = 1 channel, where the pole is located at threshold. In our derivation 
it is compensated by the same zero in the t^i channel. That is not the case of the J = 
amplitudes and therefore we are neglecting the contribution of their residue. Consequently, 
our amplitudes are not valid to obtain Adler zeros and that will affect our results at low 
energy (but no more than 0{p^)). That is another reason to differentiate the parameters 
obtained from our fit from those of the pure chiral expansion^. 



3.2.3 Multiplying by tf] 

This is apparently a harmless assumption in the above reasoning, although it dramati- 
cally affects the results of the lAM. In fact, it can happen that tf] = 0. In the (/, J) = 
(0, 0), (1, 1), (2, 0) channels of tctt scattering or in the (3/2, 0), (1/2, 0), (1/2, 1) in ttK, this 
only occurs for isolated values of s, at or below threshold. In particular, that means that the 
lAM amplitudes will have the same zeros as the lowest order Chiral amplitudes. However, 
every other partial wave vanishes at O(p^), for any s. As a consequence, the formula in 



Eq.|lJ is no longer valid. 

Nevertheless, we can generalize our previous derivation, in order to include those channels 
whose leading order is 0{p'^). We only have to go through the very same steps, although now 
we would write a dispersion relation for tjj . But let us remember that the main improvement 
of the approach is that we are calculating exactly the integral of ImG'(s) over the right cut. 
However, for that purpose we need an imaginary part, and by looking at Eq.| we can see that 
tfj = implies that Imt^j = Imt^j = 0. Therefore, unless we have a calculation up to 0{p^), 
the corresponding imaginary part will vanish. Hence when following the derivation of the 
lAM if tfj = the best we can get is plain ChPT again. At present, only 0{p^) calculations 
are available and we can only expect to obtain a real improvement with our approach in the 
six channels listed above. Thus, we will not be able to reproduce the /2(1200) resonance. 

^ While we were revising this paper it has appeared a work by M.BogUone and M.R.Pennington |Q where 
they propose other schemes with better approximations to the left cut and also include possible contributions 
from Adler zeros. 
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3.2.4 Elastic unitarity 

In order to obtain IniG on tlie rigtit cut, Eq.|T2|, we liave just made use of tlie elastic unitarity 
condition of Eq^. However, the right cut is composed of many superimposed cuts, each one 
corresponding to a different inelastic intermediate channel. Actually, Eq.^ is only true below 
the first inelastic threshold, and the real unitarity condition reads 

Imtaf3^af3 =J2^a\ ^aP^A ? ©(s - Sa) (18) 
A 

The sum is over all the physically accessible intermediate states A, whose phase space is aap- 
As far as we are neglecting electromagnetic interactions, the first inelastic channel in tttt 
is the four pion intermediate state, at 550 MeV. Similarly, for ttK is nKnir, whose threshold 
is ~ 910 MeV. Strictly speaking, only for lower energies the elastic approximation is exact. 
Nevertheless, the contribution of these intermediate states is strongly suppressed by the four 
particle phase space and we expect the lAM to provide a good approximation. 

Unfortunately within the range of energies we are interested in, there are intermediate 
channels which are not suppressed by phase space. Indeed, at approximately 985 MeV the 
inelastic KK threshold opens up. Its phase space factor is the cTq,/? in Eq.^], with Ma = Mp = 
Mk- Therefore, above the two kaon threshold we have to reconsider the derivation of the 
lAM. Let us illustrate with vrvr scattering how inelastic effects modify our result. 
As the starting point, for s > s^R) we have a new unitarity relation: 

Imt = a^^\t p +aKK I ^ (19) 

where we have denoted by t the generic tu pion elastic scattering amplitude and by tx the 
IJ partial wave of the process tttt — ^ KK. Thus we now have, for s > s^R) that 

T ^ ,(0)2 Imt/J (0)2 / , I P\ 

ImG = -t)j ^— ^ = -t\j I + aKK^-^ 1 (20) 



which differs from Eq.l2 in the term coming from two kaon intermediate production. If we 



follow the very same steps of our previous derivation, we arrive at 



tjj TT J4Af2 s'-^ys' — s — ie) \ \ t{s' 



A(s') 

Notice that, using ChPT, A(s') ^ + 0{p^). But at these high energies that is not 
negligible. Besides, we are interested in the above integral for physical values of s and 
therefore the denominator will be almost divergent for some s'. For these reasons we cannot 
neglect this integral and then we should not trust the lAM since it could miss some relevant 
physical features. 

That is indeed the case in pion scattering since, as it can be seen in Table 1, there is 
one resonance, the /o(980), whose nature is closely related to the KK threshold. Nowadays, 
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the interpretation of that resonance is still controversial: different authors propose different 



poles (not always just one) in the vicinity of the KK inelastic cut p6|, |2^. As we will see 
later, our approach is not able to reproduce any of these poles, consistent with the fact that 
the lAM makes use just of elastic unitarity. 

At this point we want to remark the importance of understanding why and when the 
method does no longer yield the right results. Let us remember that we are also thinking 
in possible applications of this unitarization procedures to the electroweak chiral effective 
lagrangian, whose reference model is the Standard Model with a heavy Higgs. In such case, 
one would expect to see a broad resonance in the scalar channel and we want to have a 
unitarization procedure whose predictions we can trust. 

3.2.5 0{p'^) approximation 

Throughout the derivation of the lAM we have been using the chiral amplitudes up to 
O(p^). Nevertheless, it is possible to extend the argument to include higher order terms, as 
for instance the 0{p^) contributions. In that case we would have started from a four times 
subtracted dispersion relation for the two-loop calculation. Once more, the integral over the 
right cut would be related to the one for G{s) = tf]'^ /tjj. Working out the expansion of the 
subtraction constants, we would then arrive to 



t 



(0)2 
IJ 



- M) ,(1) , ,(1)2 /,(o) ,(2) y^^) 

''I J '^IJ "T ^/J /''/J ''/J 

Again that is the formal [1,2] Pade approximant, and it satisfies the elastic unitarity condi- 
tion. 

As we have already mentioned two recent papers have appeared with 0{p^) calculations 
of vrvr scattering within SU{2) ChPT [|, ^. We have not used these results, since, as we 
have just seen, they will not help us to overcome any of the preceeding objections to the 
lAM. However it is quite likely that, have we used them, the parameters of the fits that we 
will present in the next sections would had been slightly modified. 



4 7T7T scattering in SUi2) ChPT 



The inverse amplitude method was first applied g, |T5| to tttt scattering without the strange 
quark. In that case, the massless limit displays spontaneous symmetry breaking from 
SU (2) L X SU{2)ji to SU{2) i^r, which is nothing but the usual isospin. The C(p^) expression 
for TTTT scattering was obtained in 0, and it is written in terms of four phenomenological 



parameters I1J2J3J4 as well as the mass and pion decay constants, and F^,. In this 
section we will review how the method is able to reproduce the p resonance. We will show 
some results for recently proposed new parameters in order to test the lAM predictive power, 
but we will also present a unitarized fit to the data. As a novelty we will use not only the 
J = phase shifts, but also those with J = 2, in order to obtain the best fit with the lAM. 
In this new calculation, we have also estimated the error bars of the unitarized parameters. 
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4.1 Results using low-energy parameters 

Let us now illustrate what happens if we apply the JAM on the ChPT amplitudes using the 
chiral parameters obtained from low energy experiments. We want to see quantitatively to 
what extent the main physical features are reproduced. 

In order to simplify the comparison with previous works, we have chosen = 139.57 
MeV and F^^ = 93.1 MeV. The values of the chiral parameters are not so clear, since they 
have considerable error bars. In Table 2 we have listed the different sets of parameters that 
we have taken from the literature to obtain Fig.l. 



Method 


h 


h 




ChPT 


-0.62±0.94 


6.28±0.48 


No resonances 


Inverse 


-0.62±0.94 


6.28±0.48 


715 MeV 


Amplitude 


-1.7±1.0 


6.1±0.5 


675 MeV 



Table 2: Sets of parameters and methods used in the text. Those in the first two fines 
come from Ku decays [ p9[ . Those in the third, from data on Kn and tttt together with some 
unitarization procedure ref. |^^. Mp is calculated with the central values. 

Let us remark at this point that for the ChPT phase shifts we are using the definition 
6 ~ (T(t(o) + Ret(i)) suggested in 0. Of course, ChPT is just a low energy approach, but 
incidentally, these phase shifts coincide with those obtained from the K-matrix unitarization 
defined as 

+ (0) I Rp+{1) 

1 - ia{m + RetW) ^ ' 

It can be easily verified that satisfies elastic unitarity, Eq.^, exactly. Consequently, the 
dotted lines in Fig.l not only give the ChPT predictions, but also the results of the K-matrix 
unitarization. We will thus confirm that such a method is not able to reproduce resonances 
by itself. They have to be added by hand. 

In Fig.l it can be clearly seen, in the J = 1, J = 1 channel, that the lAM yields a p-like 
resonance. The value of its mass is obtained from the point where 5 = 90° and it lies 10% to 
15% away from its real value. In this way, the existence of the p resonance can be regarded 
as a prediction of the lAM with ChPT and the parameters obtained from some low energy 
data. 

It is also evident that the fit of the / = 2, J = channel is correct up to much higher 
energies. In Table 2 we have also included the values of Mp corresponding to each choice of 
parameters. For all the cases we have set = 2.9 and /4 = 4.3 following reference 

The only feature of tttt scattering that is evidently missing from the unitarized results is 
the /o(980) resonance in the J = 0, J = channel. In the previous section we saw that this 
fact is connected with the failure of the whole approach to reproduce the kaon inelastic cut. 
But let us first obtain a better fit to the data. 
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Figure 1.- Phase shifts for tttt tttt. The dotted curve is plain ChPT with the U in the first 
Hne of Table 2. The other two curves are both the result of the lAM: the dashed one has been 
calculated again with the same parameters whereas the continuous one corresponds to the Ti in 
the third line of Table 2. The data come from: |l| (A), || «>, □), @ (x),g (o), H (<), 
Q (*) and (•). The results with 5L/(3) ChPT would have been exactly superimposed on 
these curves. The straight line stands aX 5 — 90°. 



4.2 Unitarized fit 

Now that we have an amphtude that describes the right cut, while keeping at the same time 
the correct polynomial form from ChPT, it seems natural to use Mp [|], |15| to fit the data. 
Note that fixing the correct mass does not imply a good fit. For instance, we could get a 
wrong width. In order to differentiate the parameters thus obtained from those coming from 
plain ChPT we will call them /i, /2. 

The (1,1) channel is almost only sensible to li — l2- With I1 — I2 = —5.95 ±0.02 we get the 
Mp listed in Table 1 and in Fig. 2 it can be seen that the results are remarkably successful. 
Later we will show that we also get the right width. 

Once that difference is fixed, we just have to determine one parameter, say l2- In previous 
studies [|, |I5| the unitarized fit to the other phase shifts was used in order to estimate the 
values of li and l2- But, as we commented above, the data in the (0,0) channel is not as 
good as that of (1,1). The same happens for the (2,0) channel, where the curves are not very 
sensible to small variations in the li parameters. Therefore, in the present work, we have also 
used the J = 2 channels (mainly that with / = 0) to further constrain the parameter range. 
Let us remember now that in these channels = and, as we have already discussed in 
Section 2, the lAM leads again to plain ChPT. That is why we will only use for them data 
up to ~ 600MeV, although in other channels we are using data at higher energies. 
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Thus, the values given in Table 3 are just a conservative estimate of the range where we 
obtain a reasonable fit in the (J, J) = (0, 0), (2, 0), (0, 2) and (2, 2) channels, when Mp fixed 
to its actual value. The results are shown in Fig.2, where the continuous line corresponds to 
the central U values and the shaded area to their uncertainties. Notice that the shaded area 
has always been obtained by varying I2 within its estimated error. 

In Fig.2 it can be seen how it is not possible to fit the /o(980) resonance with SU{2) 
ChPT and the JAM. It is clear that, even though the actual value of the 5oo phase shift may 
not lie very far from the unitarized prediction, the qualitative behavior of the curves in this 
channel is not correct above SOOMeV. 

With the li fit we can obtain the total Breit-Wigner width of the p resonance from: 

- s 

rp=^^tan5n(3) (24) 

Indeed we have computed it for different values of s around M^. The result is given in 
Table. 3 and it is quite close to the experimental result (see Table 1) although slightly high. 
We will see that it is possible to obtain the right value when using SU (3) ChPT. 

As we have already commented, this result is not at all trivial, since fitting the right mass 
does not ensure a correct description of the resonance. Therefore, even though we are now 
using the Mp experimental value, the Tp width is again a prediction of the lAM. In contrast, 
in a unitarization scheme where one introduces the resonances by hand, one has to give both 
the masses and the widths. 



Method 


h 


h 


Mp (input) 




Inverse 
Amplitude 


-0.5 ± 0.6 


5.4± 0.6 


768.8± 1.1 MeV 


155.6±1.8 MeV 



Table 3: Parameters and results of the one-loop lAM when Mp is fixed to its actual value. 



5 SU(3) Chiral Perturbation Theory 

The extension of the ChPT approach to include the strange quark was done, once more, by 
Gasser and Leutwyler . In this case there are eight Goldstone bosons, which are identified 
with the three pions, the four kaons and the eta. In principle it is possible to calculate the 
amplitudes of any process involving any combination of these particles. But the thresholds 
for these reactions are much higher than in pion scattering, which in practice restricts severely 
the effectiveness of the approach. 

Nevertheless, the lowest two particle threshold apart from tttt scattering is that of irK 
elastic scattering at 630MeV, which is still within the applicability range of ChPT. The 



calculation of this amplitude to 0{p ) was performed by Bernard, Kaiser and MeiBner ^ 
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Figure 2.- Pion elastic scattering phase shifts 5i,j obtained from the I AM fit to the correct 
Mp. The shaded areas cover the error bars of the fitted parameters with the constraint li~ h = 
—5.95 ± 0.02. The dotted straight hnes stand at S = 90°. Remember that the J = 2 partial 
waves have to be calculated as in plain ChPT. Indeed, the dashed lines in those channels 
correspond to plain ChPT with the parameters in the first row of Table 2. The symbols for 
the experimental data are the same as in Fig.l. The corresponding curves within SU{3) ChPT 
would almost superimpose. 
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who also gave the 0{p'^) result for vrvr within SU{3) ChPT. In the literature, these formulae 
have sometimes appeared with some minor errata which have been corrected in the DA$NE 
Physics Handbook pSl. However, even those formulae do not satisfy perturbative unitarity 



(see Appendix). Following the work in ||22| we have rederived an expression which does satisfy 
that requirement, and we have included it in the Appendix, together with a discussion on 
how it is obtained and its unitarity properties. 

In the SU (3) case there are more phenomenological parameters that we have set to: 

Mk = 493.65MeV , = 548.8MeV , Fk = 1.22F^ , F^ = 1.3F^ (25) 

There are also twelve one- loop parameters, denoted by L^fi). However only L^, L2, L^, L^, L 
and Lg appear in ttK in scattering, whereas in pion scattering only the following combina- 
tions are present: 

2LI + L3 LI (26) 
2LI + LI 2LI + LI (27) 

Again, and in order to simplify the comparison with previous works, we have fixed the 
following values 0]: 

LliMr,) = LliMr,) = 0.0022 L^(M^) = L^(M^) = 0.0011 (28) 

A precise value of these parameters is not very important since they are related to the 
different masses and decay constants that we had already fixed. Hence, in practice, the only 
relevant parameters for vrvr and irK scattering in SU{3) are L[,L2 and L3. 



The lAM was first applied to SU{3) ChPT by the authors in |]T3|, were we showed that it 
reproduces not only the p(770) resonance but also the i^'*(892). Our aim in this section is first 
to study the predictive power of the method, whether it can accommodate further resonant 
states, or why it cannot. Then we will present a simultaneous fit to vrvr and tcK scattering to 
the p and K* masses. The new features of this analysis is that it uses the corrected ChPT 
expressions for ttK scattering which now satisfy perturbative unitarity (see the Appendix) 
and the fact that we also use the data on the J = 2 vrvr scattering channels. We will also 
estimate the error bars on the best fit that will be used to obtain numerical values for some 
interesting phenomenological quantities. This fit will also allow us, in section 6, to perform 
a numerical study of the analytic structure of the lAM amplitudes in the complex s plane. 



5.1 Results using low energy parameters 

Let us then start with the lAM using parameters obtained from low energy data. In Table 4 
we list different choices of parameters and methods together with their results for the p and 
K* masses. As in the case of SU{2) ChPT the lAM is able to predict from low energy data 
the existence of both resonant states. Remarkably, the masses thus obtained lie again 10% 
to 15% away from their actual values. 
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Figure 3.- Phase shifts for elastic ttK scattering. The dotted curve is plain ChPT with the Li 
parameters in the first line of Table 4. The other two curves are both obtained from the lAM: 
the dashed one again with the same parameters and the continuous one with those in the third 
line of Table 4. The experimental data come from:[|8) (•), |39) (*), [|o| (o), ^ (^), ||| (□) 
and HI (A). The straight dotted line stands at (5 = 90°. 



In Fig. 3 we show the result of applying the I AM to irK scattering, with the parameters 
given in Table 4. In contrast with plain ChPT (or the K-matrix unitarization method, 

since they yield the same phase shifts), it is evident that the lAM not only accommodates 
the K* resonance, but it also reproduces the (3/2,0) channel. 

We do not display the results for vrvr scattering in SU{3) because they will almost su- 
perimpose with those in Fig.l. Indeed, the U parameters in lines 2 and 3 of Table 2 were 
obtained, respectively, from the L\,L2, in lines 2 and 3 of Table 4 ^|, by means of 



h = 9Q7r'[ALl{M^) + 2Ls- 
h = 48vr^ (4L^2(M,) - ^ - ^) 

As a matter of fact, we have calculated independently the vrvr elastic scattering in SU{2) 
and SU{3). Using the above equations to relate the parameters in both cases, and below 
kaon threshold, we have obtained the same results up to numerical differences (~1%), which 
would be unobservable in the figures. That is a nice check of our programs. Therefore, Fig.l 
is also the result for vrvr scattering in the SU{3) formalism, but now with the parameters in 
Table 4. 
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Method 


L[{Mr,) ■ 103 


L^(M^) ■ 103 


L3 ■ 103 




Mk* 


ChPT 


0.65±0.28 


1.89±0.26 


-3.06±0.92 






Inverse 
Amplitude 


0.65±0.28 


1.89±0.26 


-3.06±0.92 


717 MeV 


847 MeV 


0.6±0.3 


1.75±0.3 


-3.5±1.1 


680 MeV 


804 MeV 



Table 4: Different sets of parameters and methods used in the text. Those of the first two lines 
come from K14 decays |2^. Those of the third hne come from data on K14 and tttt together 
with some unitarization procedure (for details see ref.Q). The quoted values of Mp and Mk* 
are calculated with the central values. 



5.2 Unitarized fit 

Again we have an expression for the amphtude that behaves correctly with respect to uni- 
tarity and that presents the right form in the low energy limit. Therefore, we can try to 
use the actual p(770) and i^'*(892) masses in order to fit the rni and nK phase shifts. We 
remark once more that nothing ensures that fitting the right masses will give us the right 
description, since, among other things, the widths of the resonances could be wrong. 

When dealing with the SU{3) chiral lagrangian we have more parameters, and the way 
they appear in the amplitudes is more complicated. Let us first start with the tttt scattering 
partial waves in SU{3). As we have commented in section 3.2.1, in order to avoid confusions 
with the ChPT low-energy parameters, we will denote the parameters of our fit by L[. 

The (1,1) channel only depends on 2L[ + L3 — L2, and will be fixed with Mp. In so doing 
we get 

2L^ + L3-L^ = (-3.11 ±0.01)10"^ (30) 

As a consistency check we see that it is within a 1% of —3.14 10~^ which is obtained from 
the k parameters of the SU (2) case, with the help of Eq.p^. 

Once again we use the channels (/, J) = (0, 0), (2, 0), (0, 2) and (2,2) to determine the 
best value, which indeed is the same that we would have obtained from the I2 SU{2) 
parameter by means of Eq.p9|. It can be found in Table 5. Hence, the best SU{3) fit of the 
TTTT phase shifts, yields almost the same results as those obtained with SU{2) and the very 
same Fig.2 remains valid for 5*^7(3). Nevertheless, when computing the Tp within the SU{3) 
formalism, we obtain a much better value than in SU{2), which was about 5 MeV too high. 
It is also listed in Table 5. 

Finally we will use I/3 to fix the correct K*{892) mass. However, the i^*(892) has an 
added subtlety, namely, that the mass splitting between different charge states is of the order 
of 5 MeV. This is a small isospin breaking effect that we have not included in our approach. 
Therefore, we have used an average mass Mk* = 894.0 ± 2.5 MeV with an error bar that 
includes the mass of any K*{892) state, no matter what its charge may be. That uncertainty 
has also been taken into account in the L[ error estimates. 

Once we have L3, we use L2 and Eq.^ to obtain L^. The parameters of this fit have 



17 



Method 


L^(M^) ■ 10^ 


L^(M^) ■ 10^ 


L3 ■ 103 


Tp 




Inverse 
Amplitude 


0.41±0.20 


1.48 ± 0.33 


-2.44± 0.21 


149.9± 1.2 MeV 


41.2 ± 1.9 MeV 



Table 5: Parameters and results of the SU{3) lAM when Mp = 768.8 ± l.lMeV and M k* = 
894.00 ± 2.5 MeV are fixed to their actual values. Notice that for if* (892) we have chosen an 
average mass between its different charge states. 



been collected in Table 5, together with Tp and F^*, which can be considered as predictions 
of the approach. Notice, however, that in this case the width of the K*{892) resonance lies 
20% away from its actual value, which nevertheless is a reasonably good result in view of 
the whole fit in that channel. 

Concerning the Li parameters, they are compatible with those in Table 4, which were 
obtained from low energy data. Even more, they are also consistent with other parameters 
obtained from the lAM applied to the form factors of the K -k-kIv decays 0], which are 
very well know experimentally: 

L\{M^) = (0.74 ± 0.14)10-3 L^(M^) = (1.07 ± 0.18)10^3 Ls{M^) = (-2.45 ± 0.52)10^3 

(31) 

(notice that in that reference they are using Fx = -F^, so that the parameters do necessarily 
differ) . 

Nevertheless, it would not make any sense to try to reduce the error bars of these param- 
eters. We consider that the approach that we have been following here can only be consistent 
within a few percent error level. In order to have a better accuracy it would be necessary 
to take into account higher order ChPT corrections, isospin-breaking effects and the whole 
approach should be modified following the comments that we made in previous sections. 

In Fig. 4 we show the results of the SU (3) lAM fit to the resonance masses, in terms of 
elastic scattering phase shifts, which we think deserve some comments: 

• First notice that we are not showing the curves for tttt scattering because they are 
exactly those in Fig. 2. The differences only appear above the two kaon threshold, 
since in the SU{3) formulae we are also considering internal loops of kaons and etas. 

• In the irK —>■ ttK case we can extend the graphs up to 1100 MeV, or even more. The 
reason is that the first two body inelastic threshold is Kt] production at 1040 MeV and, 
in contrast to the tttt case, there is no nearby resonance. Indeed, the next resonant 
state in irK elastic scattering is ^^0(1430), very high to affect dramatically our results 
at 1100 MeV, but also to be correctly reproduced by the lAM method. Nevertheless, 
the existence of the Kt] threshold can be noticed in the I = 1/2, J = channel, as a 
small bump in the curves at precisely 1040 MeV. 
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Figure 4.- ttK elastic scattering phase shifts Sij obtained from the lAM fit to the correct Mp 
and Mk'- The shaded areas cover the error bars of the fitted parameters with the constraint 
2L\ + L3 - = (-3.11 ± 0.01) 10"^ The dotted straight hne stands ai 5 ^ 90°. The symbols 
for the experimental data are the same as in Fig. 3. 



• The shaded area in the i^*(892) channel is not only due to the averaged mass for 
i^*(892) with 2.5MeV error, but also to the fact that we have to determine several 
parameters to get the right mass, in contrast with the p(770) case, when we only had 
to fix one. 

We have explicitly checked that our ChPT amplitudes satisfy perturbative unitarity. 
As it is explained in the Appendix, previous calculations [^, [1^, including ours, did not 
respect this condition, although by a very small amount. That is why the values of the best 
parameters for this fit are slightly different from those of our previous work [O . 



Phenomenological parameters 



Once we have a good parametrization of tttt and txK elastic amplitudes, we can use it to 
obtain the values of some relevant phenomenological parameters. First we can calculate the 
scattering lengths, which determine the strength of the interactions at low energy. Despite 
our lAM fit makes use of high energy data, we expect that it will reproduce the low energy 
behavior since in the low energy limit it reduces to the chiral expansion, which at 0{p'^) 
already yields quite good values (see Tables 6 and 7). However, as far as the JAM is non- 
perturbative we are also taking into account higher order effects, that will modify the results. 
Indeed, some of these lengths have already been calculated with the lAM and it yields slightly 
better results than plain ChPT [Q. We have made again the calculation with our fit, but 
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as far as we have an estimate of the error bars in the Li parameters, we will also give the 
error estimates coming only from the uncertainties in Lj (mostly dominated by that of L2). 

Before giving the results, it is convenient to recall that the scattering lengths have two 
different normalizations. Namely: 

^etjj{s) = q^' (a'j + + 0{q'')) (32) 
for TTTT scattering, where q is the CM. momentum q^ = s/A — M^, and 

RetrAs) = ^q'' (a^ + + Oiq')) (33) 

for ttK scattering, where now q"^ = [s - {Mk + M^y][s - {Mk - M^'f]/4s. 

The predictions of our fit for the vrvr and irK scattering lengths are given in Tables 6 and 
7 (in units). Notice that all the values are compatible with the experimental data, and 
in general they only differ very slightly from the 0{p'^) ChPT results, usually in the right 
direction toward the central value. However, the experimental error bars are still too big to 
arrive at any conclusion. Also the error bars in the lAM have to be interpreted cautiously, 
since they are obtained only from the uncertainties in the Lj parameters. 

As we have already commented, very recently there has appeared a two loop calculation 
of TTTT scattering within SU(2) ChPT. It estimates a[j ~ 0.217 or 0.215 and a[j -a^ ~ 0.258 or 
0.256, which are precisely the values obtained with our lAM fit. This fact gives support to 
the idea that the lAM somehow takes into account higher order terms even at low energies. 

Notice that we do not compare with the two-loop calculation of the scattering lengths and 
slopes in because they have used them as an input in a x^-fit to determine their additional 
a and f3 parameters. Therefore, their values are almost exactly those of the experimental 
data. However, as far as there is no data for the bl im slope parameter, their value can be 
regarded as a prediction. They give b\ = (0.54 ± 0.15)10"^, which is consistent with our 
result and with b\ = (0.6 ± 0.4) lO'^, that was obtained from sum rules in [Q. 



We have also calculated the phase of the e' parameter, which measures direct CP violation 
in K —>■ TTTi decays [^. It is related to the s-wave phase shifts as follows: 

0(e') = 90° - (5° - 6l),=,p^ (34) 

Our result is: 

m = (42;D° (35) 

very close to <^(e') = (45 ± 6)° which is obtained in plain ChPT |^^. In contrast with the 
case of the scattering lengths, the value of this angle is not used as an input in Q and is 
therefore a prediction of their best fit. The value they quote is 0(e') = (43.5 ± 2 ± 6)". 

Finally, in Fig. 5 we show the phase difference 600 — Su, compared with the hitherto 
available experimental data [^. The difference between the lAM and plain ChPT at high 
energies is due to the presence of the p resonance. Nevertheless, there are also some differ- 
ences at low energies, since the dispersive approach is somehow taking into account higher 
order contributions. 
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ChPT 


lAM fit 


Experiment 


at 


0.201 
0.26 


0.216 ± 0.008 
0.289 ± 0.025 


0.26 ± 0.05 
0.25 ± 0.03 


al 
bl 


-0.041 
-0.070 


-0.0417 ± 0.0014 
-0.075 ± 0.003 


-0.028 ± 0.012 
-0.082 ± 0.008 


a\ 
bl 


3.6 • 10-2 
0.43-10-2 


(3.744 ± 0.002)-10-2 
(0.515 ± 0.001)-10-2 


(3.8 ± 0.2)-10-2 




20 ■ 10-^ 


(17.1 ± 3.5)-10-^ 


(17 ± 3)-10-^ 


al 


3.5 ■ 10-^ 


(2.8 ± 1.5)-10-^ 


(1.3 ± 3.1)-10-^ 



Table 6: rnr scattering lengths. The one-loop ChPT results are taken from |2^. The experi- 
mental data come from . The errors in the lAM fit come only from the uncertainties in the 
parameters . They do not include other theoretical uncertainties. 



a^j ChPT lAM fit Experiment 



3/2 


-0.043 


-0.049 ± 0.004 
-0.026 ± 0.003 


-0.13.. 


.-0.05 


1/2 
% 


0.148 


0.155 ± 0.012 
0.087 ± 0.016 


0.13.. 


.0.24 


1/2 

ai 


0.012 


0.0146 ± 0.0012 


0.017.. 


.0.018 



Table 7: ttK scattering lengths. Note that the ChPT results have been obtained using the 
corrected formulae in the Appendix. The experimental data come from |22] 
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Figure 5.- 600 — 611 phase shift difference from the lAM fit (sohd hne) and plain ChPT (dashed 
hne). The shaded area covers the uncertainty in the parameters and the data come from 

6 The lAM in the complex s plane 

The main objection to unitarization procedures is the apparent arbitrariness in their pre- 
dictions, which may differ from one another. In most cases, these methods are nothing but 
a small modification of the amplitudes so that they can satisfy the unitarity constraint in 
Eq.^, while keeping at the same time the good low energy behavior. But that constraint 
is not enough to determine the amplitude completely. Thus there are as many unitariza- 
tion techniques as algebraic tricks to implement such a constraint exactly or to get a better 
approximation. 

However, we have already seen in Section 3.1 that, below any other inelastic threshold, 
the inverse amplitude method can be derived directly from the analytic structure of the 
general two body elastic scattering amplitude. Our purpose in this section is to show that, 
apart from satisfying elastic unitarity, it provides the correct analytic structure required 
from relativistic Quantum Field Theory. Such an structure is not trivial at all and cannot 
be reproduced by other unitarization procedures. Both the left and right unitarity cuts are 
already present in plain ChPT, therefore, we will mainly focus on the poles in the second 
Riemann sheet. 

In the previous section we used the most naive criteria to identify resonances, i.e., that 
the phase crosses the 6 = 90° value. However, that is only true for the simplest cases. The 
rigorous characterization of resonances is made in terms of poles in the second Riemann 
sheet of the amplitudes in the s complex plane. Indeed, when a resonance is produced by 



22 



just one of these poles, both its mass and width are related to the pole position by 



Sp,,,^MR + ^^^ (36) 



provided the width is small enough. 

In this work we have extended to the s complex plane both the vrvr and the nK elastic 
scattering JAM amplitudes obtained in the previous section. Notice that the cuts in ChPT 
come from logarithmic functions, so that we have infinite sheets in the complex plane. How- 
ever, only two of them correspond to the first and second Riemann sheets. Once we have 
identified these sheets we can check whether the resonances that we found in previous sec- 
tions are produced by a pole in the second Riemann sheet and thus whether they have a real 
sound basis. 

We will first analyze the tttt hit process. In Figure 6 we represent the imaginary part 
of the amplitude in the complex s plane for the three channels [I, J) = (0,0), (1,1) and 
(2,0). Notice that when we say complex s plane, we mean that we have parametrized s as 
s = {E + iCy, where E is the CM energy and is represented in the real axis whereas C 
provides the complex part. On the left column we have displayed the results in the first 
Riemann sheet, whereas in the right column we have continued through the cut to the lower 
half of the second Riemann sheet. In all cases it can be clearly noticed the existence of a 
cut on the real axis on the first Riemann sheet. As we had commented before, a right cut is 
not anything completely new, since it is already present in one-loop ChPT, although in that 
case, the values that the amplitudes take on it are different. In contrast, the most striking 
new feature in the lAM amplitudes is the appearance of poles in the second Riemann sheet 
and how they determine the amplitude shape for the physical values of s. 

Indeed, we have found two poles with Ims < in the second Riemann sheet, one in the 
(0,0) partial wave and another one in (1,1). Let us start with the second, which clearly 
corresponds to the p resonance. The position of this pole can be obtained from the contour 
plots in Fig. 7, and it is found at around Eji ~ 760 — i75. Using Eq.^ we see that it is in 
a good agreement with the p(770) mass and width parameters given in Table 1. Therefore, 
we can conclude that this pole is completely consistent with the p(770) resonance. 

The other pole that can be seen in vrvr — tttt is on the (0, 0) channel. Using the parameters 
of the best SU{3) lAM fit of the previous section, we find that it is located at Ef> ~ 440— z245. 
It is not responsible for the appearance of any resonance, since it is very far away from the 
real axis. However, from purely phenomenological fits to pion scattering data it had already 



been pointed out the existence of such a pole around Eji ~ 408 — i342 MeV [0. This pole 
is responsible for the strong interaction in that dominates the at low energy the two pion 
(0, 0) channel. We can now see that even in the channel where there is not an apparent 
improvement, the lAM yields the correct analytic structure. 

Much as it happened in previous sections, the method is not able to reproduce the /o(980) 
resonance. As we already commented, the interpretation in terms of poles of this resonance 
is still controversial. Following the same steps as before, we have also identified the four 
Riemann sheets that now appear due to the superposition of two cuts. Indeed, we have even 
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Figure 6.- Imaginary parts of the tttt — > tttt amplitudes in the complex s plane. The first 

row is the (/, J) = (0,0) chamicl, the second is (1, 1) and the bottom is (2,0). The left plots 
correspond to the first Riemann sheet, and those on the right, to the second. 



24 



700 710 720 730 740 750 760 770 780 



380 390 400 410 420 430 440 450 



850 850 870 



890 900 910 920 930 940 950 



Figure 7.- Contour plots of the second Riemann sheets for different SU{3) ChPT unitarized 
amphtudes. From left to right they correspond to the (1, 1) and (0,0) tttt scattering channels 
and the (1/2, 1) ttK ttK channel. 



implemented the lAM derived with the inelastic unitarity condition in section 3.2.4. We 
have not found any pole that could hint at the existence of such a resonance. 

As we have already explained, we should not expect to find anything since the approach 
is not able to reproduce properly the two kaon unitarity cut and consequently neither its 
associated sheet structure. 

Let us now address to the irK elastic scattering case. Again, in Fig. 8 we have displayed 
the imaginary part of the amplitudes for the (3/2,0), (1/2,0) and (1/2, 1). Those pictures 
on the left represent the first Riemann sheet and those on the right, the second. Once more 
the existence of a unitarity cut is clear, but there is also the appearance of a pole in the 
appropriate channel. In particular, using the third contour plot of Fig. 7, we have found a 
pole in Eji ~ 890 — 'j20 MeV, which using Eq.^ yields again the mass and width for the 
i^*(892) resonance that we gave in Table 5. 



7 Conclusions 

In this work we have shown how the lAM provides a consistent technique to accommodate 
resonances. Indeed, based on its derivation from Dispersion Theory, we have made a system- 
atic analysis of its applicability, which is mainly limited by the existence of two body inelastic 
thresholds and by the fact that the tree level approximation vanishes in some channels. 

We have found that it is able to predict the most relevant features of strong elastic 
scattering, once the chiral parameters are determined from low energy data. Quantitatively, 
the errors are hard to estimate, but we have found in all cases that the mass of the predicted 
resonances fall within approximately 15% of their actual values. We think this fact gives 
a sound basis for its application in order to obtain at least a qualitative description of 
resonances in the strongly interacting symmetry breaking sector. 

Moreover, once we force the lAM results to fit the actual resonance mass values, we get a 
remarkably good fit which is able to reproduce the experimental data up to the next relevant 
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Figure 8.- Imaginary parts of the -kK — > -kK IAM amplitudes in the complex s plane. The first 

row is the (/, J) = (3/2,0) channel, the second is (1/2,0) and the bottom is (1/2,1). Again, 
the left plots correspond to the first Riemann sheet, and those on the right, to the second. 
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two body inelastic threshold. Following that procedure, we have given the unitarized SU{2) 
ChPT fit to TTTT — > 7171 as well as that of SU(3) to tttt and ttK elastic scattering. For the 
first time we have estimated the values of the unitarized chiral parameters together with 
their error bars. These values do not lie very far from those obtained without the lAM, 
and therefore do not spoil the low energy expansion, as can be noticed from the scattering 
lengths that we have given. 

With this fit, we have calculated several low-energy phenomenological parameters, like the 
scattering lengths. Our values differ from those obtained at 0{p'^) due to the unitarization. 
However we expect that they include other corrections due to unitarity and resonant effects. 

We consider that it would not make any sense to try to reduce the error bars in the 
unitarized parameters within this approach. One has to keep in mind that we have neglected 
higher order ChPT corrections, isospin-breaking contributions, and that we have used high 
energy data which is very sensible to such effects. It is quite likely that, in order to obtain 
results consistent to a higher degree of accuracy, the lAM in the simple version that has 
been used here, will not be enough. 

Finally we have also shown how the lAM yields the proper analytic structure in the 
complex s plane, in contrast with other unitarization techniques. Indeed, we have found 
that the apparent resonant behavior that is observed on phase shifts, is produced by the 
corresponding poles in the second Riemann sheet, meeting the strict requirements imposed 
by general relativistic Quantum Field Theory. 

Therefore, we think that the lAM and unitarization by means of dispersion theory is 
the most natural and economic way to extend the applicability of Chiral Lagrangians. We 
have seen however, that its main limitations come from the existence of two body inelastic 
thresholds. Nevertheless, work is still in progress in the subject, the lAM has been recently 
applied to other processes and higher order ChPT calculations will be soon available. As far 
as some other physically relevant features do not lie very far from the present applicability 
limits, it seems very likely that they can be reproduced in the near future. 
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A Elastic scattering amplitudes in SU{3) ChPT and 
Unitarity 

The first calculation of elastic ttK scattering was later performed in [^]. These amplitudes 
were given in terms of physical as well as lowest order masses and decay constants, which 
are usually denoted by Mp,Fp and Mp,F°, respectively (P being either 7i,K or rj). Of 
course, the only measurable parameters are the first, and when comparing with experimental 
observations, one has to eliminate those from lowest order in terms of the physical constants. 
Indeed, it is possible to find the relation between M° and M^r as well as that between 
and M^- Unfortunately, at lowest order there is only one which is related both to 
and Fk- Hence, whenever one finds F° in an expression there are two choices: either 
relate it to F.,^ or to Fk- The difference between the two choices will be one order higher in 
the chiral expansion. For instance, if one has an O(p^) expression with Fq, in principle one 
can substitute it by F^, F^ or ■ F^- AH these choices are equally acceptable. When one 
is working only with pions, the natural choice is the one that leaves all the expressions in 
terms of F^^. When one is dealing both with pions and kaons, it is not so obvious. However, 
once one choice is done for the 0{p'^) term, we have to keep it for the 0{p'^) contribution, 
otherwise one would violate perturbative unitarity, Eq.^. 
Surprisingly, in the amplitude in the literature 

the choice for the 



in our previous work |]T3| 



which is the one we had also followed 
term is different from that of the 
contribution that yields the imaginary part. Indeed, the 0{p'^) term is written just in terms of 
whereas is written in terms of F^^ and F^- As a consequence, there is a F^/F^ factor 
of difference between Imt^^^ and o'^^l't^'^^l'^- Numerically that amounts to a (1.22)^ ~ 1.5 
factor. 

Thus, we have rederived from the original work p2| the amplitudes in terms of physical 
quantities, (also correcting some small errata) so that they satisfy perturbative unitarity. 
We have chosen to write the formulae symmetrically with respect to and Fk- But the 
other choices are equally acceptable. The result is: 



T'/^is,t, 



u 



Tj{s,t,u) 



T[{s,t,u) 



T^is,t,u) 



Ml + Ml-s 



+ T/(s, t, u) + (s, t, u) + Tf (s, t, u) + 0{s 



'^Ft.Fk 



(37) 



IGF^Fk 

{AL[it - 2Ml){t - 2Mi) + 2LI 



Fm 



K 



Ml - Mlf + {u-Ml 



M 



2 \2 
K) 



+ LI 



u-M^- Mf^Y + {t- 2M^){t - 2Mj^] 



t{Ml + M|) - AMlMl 



+ 2LlMl{Ml-Ml-s)+S{2Ll + Ll)MlMl] 



AF^Fl 1 2 
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[s - t) (l^k{u) + LKr,{u) - u {mIj,{u) + Ml,^{u) 



+ (Ml - mI)\m:m + M\ 
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+ ^(Mi - MI)[K^k{u){^u - 2Ml - 2Ml) + Kk^{u){?,u - 2M|. - 2Ml)] 
+ - 12«(M|. + Ml) + 4(Mi + M2)2] + +Jl^{s){s -Ml- Mf 



+ -rK,iu)(u--iMl + Ml)] +-Jl^{t)t{2t-Ml) + -J],^{t)t' 



+ \r,,ml{t-lMl)] 



The functions MpQ, LpQ, Kpq, Jpq, ^p, with P,Q = TT,K,ri, can be found in although 
they should be written in terms of physical quantities. 

We have verified analytically that this amplitude satisfies the perturbative unitarity con- 
straint. Moreover, we have used that constraint as a check of our programs. 

We want to remark again that this way to write the 7iK amplitude is one of several 
possible choices, since we could have chosen to write everything just in terms of F^^, for 
example. The important point is to keep the same choice both for the and the 0{p^). 

For completeness, we will also give the SU{3) formulae used in this work for tttt scattering, 
because they have also appeared with some minor errata in the literature: 



A{s,t, 



u 



B{s, t, u) 



l£_Jfil + B{s, t, u) + Cis, t, u) + Ois^) 
1 f M"^ 1 1 



(3^ 



+ \{t- 2Ml)'JUt) + tis - u) \Ml^{t) + ^M],^(t)] + (t ^ «)} 
C{s,t,u) = -lr,{{2L\ + L,){s-2Mlf + Ll[{t-2Ml)' + {u-2Mlf] + 



+ 



{ALl + 2Ll)Ml{s - 2Ml) + (8L^ + ALl)Ml} 
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